Onset of dendritic flux avalanches in superconducting fllms 
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We report a detailed comparison of experimental data and theoretical predictions for the dendritic 
flux instability, believed to be a generic behavior of type-II superconducting films. It is shown that a 
thermo-magnetic model published very recently [Phys. Rev. B 73, 014512 (2006)] gives an excellent 
quantitative description of key features like the instability onset (first dendrite appearance) magnetic 
field, and how the onset field depends on both temperature and sample size. The measurements 
were made using magneto-optical imaging on a series of different strip-shaped samples of MgB2. 
Excellent agreement is also obtained by reanalyzing data previously published for Nb. 



Phenomena that create intriguing traces of activ- 
ity that can be observed by direct visual methods are 
among the most fascinating things in nature. Penetra- 
tion of magnetic flux in type-II superconductors seen by 
magneto-optical (MO) imaging is one example, where es- 
pecially the spectacular dendritic flux patterns occur- 
ring in superconducting films are currently attracting 
much attention. The phenomenon has been observed in 
a large number of materials; YBa2Cu30aj, Nb, MgB2, 
NbaSn, NbN, YNisBsC and Pb [1-7], all films, and 
showing essentially the same characteristic behavior. In 
abrupt bursts the film becomes invaded by flux in narrow 
finger-like regions that often form a complex and sample- 
spanning dendritic structure. These sudden events occur 
typically during a slow ramping of the applied magnetic 
field, and at temperatures below a certain fraction of the 
superconducting transition temperature, Tc. It is also 
characteristic that the flux patterns are never reproduced 
when experiments are repeated, see Fig. 1, thus ruling out 
possible explanations based on material defects guiding 
the flux motion. The massive experimental data existing 
today [1-15] indeed suggest that the phenomenon is a 
generic instability of the vortex matter in superconduct- 
ing films. 

Abrupt flux avalanches are known to occur in super- 
conductors for two fundamental reasons: (i) the motion 
of vortices releases energy, and hence increases the lo- 
cal temperature, and (ii) the temperature rise reduces 
fiux pinning, and facilitates further vortex motion. This 
makes up a positive feedback loop that may lead to 
an instability [16, 17]. However, why such avalanches 
should develop into dendritic patterns is a topic under 
vivid discussion, and several competing theories were re- 
cently proposed. They include a stability analysis taking 
into account the complicating non-local electrodynamics 
of thin film superconductors [18, 19], a boundary layer 
model assuming shape-preserving fronts [20] , and a shock 
wave approach [21], all leading to substantially different 
predictions. In this work we report on the first experi- 




FIG. 1: (a). Three MO images of flux penetration in MgB2 
taken during repeated identical experiments, (b). Image ob- 
tained by adding the 3 complementary colored images above. 
In the sum image the grey tone regions are those of repeated 
behavior, whereas colors show where there is no or only par- 
tial overlap. Strong irreproducibility is seen in the dendrite 
shapes, while the penetration near the edge and along static 
defects is reproducible. The dendrites tend to nucleate at pre- 
ferred sites along the edge, which is due to small edge cavities 
giving local field amplification. The experiments were per- 
formed after cooling to 9.2 K and applying a magnetic field 
of 20 mT. 



ments designed specifically to check the validity of these 
models. It is shown that the model in Ref. 19 provides 
an excellent quantitative description of key features, such 
as the instability threshold field, i^th, i-e., the magnetic 
field when the first avalanche occurs, and how i^th de- 
pends on both temperature and the sample size. The 
results were obtained by MO imaging of flux penetration 
in MgB2 films. Also earlier observations of Nb films [10] 
are shown to be in full agreement with this model. 

Thin films of MgB2 were fabricated by a two step 
process [22], where first a film of amorphous boron was 
deposited on an AI2O3 (1102) substrate using a pulsed 
laser. The B film and high-purity Mg were then put into 
a Nb tube, which was sealed in a high purity Ar at mo- 
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FIG. 2: MO image showing flux distribution in MgB2 strip- 
shaped samples at 4 K and 15 mT applied field. The image 
brightness represents the local flux density. Both the number 
and size of the dendrites are larger for the wider samples. 

sphere and post-annealed at 900° C. To eliminate possible 
contamination with oxygen, water, and carbon, the sam- 
ples were not exposed to air until the final form of the 
film was produced. The MgB2 films possess c-axis ori- 
entation, as confirmed by scanning electron microscopy, 
and magnetization data show a sharp superconducting 
transition at 39 K. The film thickness was 300 nm. 

A set of eight MgB2 film samples was shaped by photo- 
lithography into 3 mm long rectangles having different 
widths ranging from 0.2 mm to 1.6 mm. All the samples 
were made from the same mother film, allowing simul- 
taneous and comparative space-resolved magnetic obser- 
vation [23]. An additional 5 mm wide sample was made 
using the same preparation conditions. A standard MO 
imaging setup with crossed polarizers and a ferrite garnet 
indicator was used to visualize flux distributions. Shown 
in Fig. 2 is an image of the flux penetration pattern when 
the eight samples, initially zero-field-cooled to 4 K, were 
exposed to a perpendicular applied magnetic field slowly 
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FIG. 3: Threshold magnetic field for onset of the dendritic 
instability in MgB2 strips of different width (symbols) plot- 
ted together with a fitted theoretical curve (full line), which 
diverges at a finite w indicated by the dashed asymptote. 

ramped to 15 mT. The magnetic flux enters the super- 
conductor in a form very much dominated by abrupt den- 
dritic avalanches, although quite differently for the vari- 
ous samples. It is evident that the number of dendrites, 
their size and branching habit depend strongly on the 
sample width. Whereas the wide strips become densely 
filled with fiux dendrites, the more narrow samples con- 
tain fewer, until at the 0.2 mm wide strip fiux dendrites 
almost never appear. 

This qualitative result, was followed up by measur- 
ing how the instability threshold field i^th depends on 
the strip width. Results obtained for all eight strips are 
shown in Fig. 3, where each data point represents an 
average over 4 repeated experiments using identical ex- 
ternal conditions. The error bars indicate the scatter in 
the observed i^th- A variation as much as 30% implies 
that the nucleation of this instability is strongly affected 
by random processes, which is also consistent with ear- 
lier experiments [2, 3, 5, 6, 8, 10]. Nevertheless, the data 
in Fig. 3 show a clear increase in the threshold field as 
the strip becomes narrower. In other words, reducing the 
sample width increases the stability of the superconduc- 
tor. 

Measurements of the temperature dependence of i^th 
are shown in Fig. 4. One sees that i^th not only increases 
with temperature, but appears to diverge at a certain 
temperature. Above this threshold temperature, Tth, 
found to be close to 10 K for MgB2 films, the dendritic 
instability disappears entirely. Included in the figure are 
also data we have extracted from a previous MO inves- 
tigation of dendritic fiux penetration in Nb films [10]. 
The two behaviors show remarkable similarities, although 
with different threshold temperatures, approximately 6 K 
in the Nb case. 

To explain these observations we adopt the model de- 
veloped in Refs. 18, 19. There, a linear analysis of the 
thermo-magnetic instability in a long and thin super- 
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FIG. 4: Temperature dependence of the threshold magnetic 
field. Experimental data obtained for our largest MgB2 sam- 
ple and for a 1.8 mm wide Nb film [10] are plotted as • and 
A, respectively. The full lines are theoretical fits. The dashed 
lines show the limiting temperature above which the instabil- 
ity vanishes. 

conducting strip thermally coupled to the substrate was 
worked out. The analysis considered a strip of width 
2w and thickness d <^ w placed in an increasing trans- 
verse magnetic field leading to a quasi-critical state in 
the flux penetrated region near the edges. By solving 
the Maxwell and the thermal diffusion equations, it was 
shown that for small fields there are no solutions for per- 
turbations growing in time, implying a stable situation. 
As the field increases the distribution can become unsta- 
ble, with a fastest growing perturbation having a non- 
zero wave vector along the film edge. This means an 
instability will develop in the form of narrow fingers per- 
pendicular to the edge - a scenario closely resembling the 
observed dendritic flux behavior. Within this model, the 
threshold flux penetration depth, £*, when the supercon- 
ducting strip first becomes unstable, is given by Eq. (25) 
of ref. 19, which can be expressed as 



Here j'^ is the temperature derivative of the critical cur- 
rent density, n is the thermal conductivity, and is the 
coefficient of heat transfer from the superconducting film 
to the substrate. The parameter n ^ 1 characterizes the 
strongly nonlinear current- voltage curve of the supercon- 
ductor, described by the commonly used relation for the 
electrical field, E oc j^. 

The threshold field, i^th, is obtained by combining 
Eq. (1) with the Bean model expression for the flux pen- 
etration depth of a long thin strip in a perpendicular 
applied field [24, 25], 

Hth = — arccosh ( — ) . (2) 



Plotted in Fig. 3 as a solid line is this function using 
jf'c = 9 X 10^° A/m^, a value obtained for MgB2 at 4 K by 
extrapolation of jc(T)-curves measured under the stable 
conditions above Tth- The only adjustable parameter, 
£*, was chosen equal to 82 /im, which gives an excel- 
lent agreement with our data. It follows from Eq. (2) 
that narrower strips need a larger field to reach the crit- 
ical penetration depth which is exactly what we find 
experimentally. Furthermore, the model predicts that 
i^th should diverge when the strip halfwidth decreases 
towards w = ^ also this fully consistent with our MO 
observations. 

To fit the observed Hi\^{T) one needs temperature de- 
pendent model parameters. We assume then a cubic de- 
pendence of the thermal conductivity, = k{T/Tc)^^ as 
suggested by low-temperature data for MgB2 [26] . Simi- 
larly, a cubic dependence of the heat transfer coefficient, 
^0 = ^0 {T /Tcf is chosen in accordance with the acous- 
tic mismatch model confirmed experimentally for many 
solid-solid interfaces [27]. Furthermore, we assume a lin- 
ear temperature dependence for the critical current den- 
sity, jc = jco(l — T/Tc), and with a pinning potential, 
U oc 1 — T/Tc, the exponent n ^ U/kT also becomes 
T-dependent, n = h {Tc/T — 1). 

Combining all these equations, one obtains a theoret- 
ical i7th(T), and shown in Fig. 4 are such curves fitted 
to both sets of experimental data. The model clearly 
reproduces the two key features; (i) the existence of a 
threshold temperature Tth above which the instability is 
absent, and (ii) a steep increase of the threshold field i^th 
when T approaches Tth- For MgB2 the fit was made with 
jco = 10^^ A/m^ and k = 160 W/Km [26], and choos- 
ing n = 10 corresponding at T = 10 K to the commonly 
used n = 30. The remaining parameters are the electric 
field and the heat transfer coefficient, where best fit was 
obtained with £; = 30 mV/m and ho = 17 kW/Km^. 
It should be emphasized that the experimental data for 
both Htii{w) and i^th(^) were fitted using the same pa- 
rameter values, and in both cases giving excellent quanti- 
tative agreement. Figure 4 also shows a similar fit for the 
data obtained for Nb, using = 9.2 K, jco = 10^^ A/m^, 
w = 0.9 mm, d = 0.5 /im, [10] k = 120 W/Km, [28] 
n = 40, £; = 200 mV/m and ho = 36 kW/Km^. Again 
the model excellently describes the experimental behav- 
ior. 

The fitted electric fields represent upper limiting val- 
ues, since we used bulk values for thermal conductivity, 
which in general are larger than for films. Nevertheless, 
both values largely exceed the estimate, E ~ H£^^ ex- 
pected for a uniform and gradual flux penetration with 
a ramp rate of H ^ 1 mT/s as used in the experiments. 
We believe this discrepancy is due to the fact that lo- 
cal, rather than global, conditions govern the onset of 
the instability. Assuming that the flux dendrites are nu- 
cleated by abrupt microscopic avalanches of vortices [29], 
local short-lived electric fields can easily reach those high 
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FIG. 5: Theoretical stability diagram predicting the threshold 
temperature Tth for different film width. The curve is plotted 
for parameter values corresponding to MgB2 films. 

values. In fact, such avalanches consisting of 10^ — 10^ 
vortices occurring in an area of ~ 20 /im were recently ob- 
served by high-resolution MO imaging in MgB2 films [30]. 
Electric fields close to 30 mV/m would be created if such 
avalanches occur during a time span of the order of 10~^ 
seconds. Randomness in such avalanches may also ex- 
plain the large scatter of the observed JYth values. We 
also note that the estimated electric field at the nucle- 
ation stage is still much lower than E values at the tip 
of an already propagating dendrite [31]. This fact is in 
agreement with expectations. 

Finally, we emphasize that the two functions i^th('^) 
and Hth{T) have a similar feature, namely a divergence 
at some value of the argument beyond which the system 
becomes stable, see Fig. 3 and Fig. 4. These stability 
thresholds are actually related to each other by the 
condition ^*(Tth) = w. The relation between the 
threshold temperature and the strip width is shown 
in Fig. 5, and represents the stability diagram in w — T 
coordinates, here plotted for parameters valid for MgB2. 
It follows from the model that the temperature range 
of the instability increases monotonously with the strip 
width, but is limited upwards by a temperature close to 
10 K for large-size films, as confirmed by many previous 
experiments [3, 8, 12-14]. The general result that the 
instability is suppressed for sufficiently narrow strips, is 
of particular importance for design of superconducting 
electronic devices or other applications making use of 
thin film superconductors operating at temperatures 
below the instability threshold value. 
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